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THE PHOTON GAS FORMULATION OF THERMAL RADIATION 


By Robert C . Ried, Jr. 
Lyndon B. Johnson Space Center 

SUMMARY 


Thermal radiation transport is formulated on the basis of a photon gas and 
also , independently , on the basis of the corresponding transverse electromagnetic 
wave formulation . Nonequilibrium statistical thermodynamics is employed in the 
form of a postulated entropy extremum principle to truncate the photon gas moment 
equations by way of the most probable distribution function subject to appropriate 
physical constraints . This approach is applied to an equilibrium statistical thermo- 
statics derivation of the energy , linear momentum , and intrinsic angular momentum 
equations for an isotropic photon gas . Conservation equations are obtained in 
terms of a temperature, an energy average velocity, and an intrinsic angular mo- 
mentum diffusion velocity . Interactions between the radiation and physical media 
are considered only in a general manner. The relativistic variation of a photon 
gas temperature is obtained , and the entropy is established as a true scalar . 

Major appendixes are included: (1) to illustrate the general form of photon- 

media interactions , (2) to provide an independent transverse electromagnetic wave 
formulation of radiation, (3) to establish the photon gas phase space representation 
by direct quantization of the electromagnetic variables , (4) to summarize integrals 
required in the statistical formulation, (5) to provide the correspondence between 
the thermodynamic constraint of momentum and a special relativistic Lorentz trans- 
formation, and (6) to illustrate the application of nonequilibrium statistical thermo- 
dynamics by the treatment of a nonisotropic photon gas with the subsequent devel- 
opment of a corresponding second-order tensor temperature . 


INTRODUCTION 


The quantitative description of thermal radiation transport is generally 
a complex problem inherently involving both the microscopic and macroscopic 
characteristics of a physical medium. At the microscopic level, classical descrip- 
tions of the physical medium and of the electromagnetic field generally are not 
adequate for the description of mutual interactions . A quantum mechanical treat- 
ment of the particles forming the physical medium and a quantum electrodynamic 
treatment of the transverse electromagnetic waves are required to describe the mu- 
tual interaction . These treatments culminate in transition probabilities for emission 
or absorption by the physical particles of electromagnetic angular momentum 
quanta having discrete energy and linear momentum , commonly called photons . 

At the macroscopic level, the radiation emission and absorption characteristics 



of the medium involve a statistical consideration of the energy , the linear momentum , 
and the angular momentum of the physical particles forming the medium . The 
objective of this report is to provide a comparable statistical consideration of 
the energy, the linear momentum, and the angular momentum of the photons 
that make up the thermal radiation field and thereby to provide a more comprehen- 
sive basis for the treatment of macroscopic radiation transport . 

Thermal radiation transport in astrophysics , physics , and engineering 
problems has generally been treated by a semimicroscopic approach through 
the direct application of a radiation transfer equation (refs. 1 and 2) . This 
approach amounts to photon counting without statistical consideration of the photons 
themselves. The radiation transfer equation can be solved, in general, along 
any particular direction in space . However , the coupling of these solutions 
for the infinite number of directions at every point in space is limited to special 
cases . In an optically thin medium , the coupling of radiation in different directions 
is relatively weak , and the direct application of the radiation transfer equation 
provides a reasonable approximation . In one-dimensional radiation transport 
problems , symmetry implicitly provides a statistical consideration of photon 
behavior , and the direct application of the radiation transfer equation is possible . 

In optically thick media , a diffusion approximation (refs . 3 and 4) to the radiation 
transfer provides a realistic approach to radiation transport within the media; 
however , this approximation generally breaks down at physical boundaries . 

An alternate approach toward obtaining a description of thermal radiation 
transport is through a moment method in which the Boltzmann equation for photons 
is expanded in a series of moments of the photon distribution function. Use of 
this method produces an infinite set of equations in physical space and time that 
must be truncated for practical application. In practice, these moment equations 
are normally truncated somewhat arbitrarily for mathematical convenience . 

This report includes a moment expansion obtained by analogy with molecular 
gas dynamics (refs. 5 to 13); hence, the photon gas designation. The moment 
equations are truncated physically by assuming that the general thermodynamic 
entropy will be a maximum subject to the constraints on the system, an assumption 
that , in principle , can be applied at any moment level . 

In general , this is a nonequilibrium statistical thermodynamics approach 
toward a macroscopic description of thermal radiation transport . The approach 
is developed generally and then applied to the restricted equilibrium statistical 
thermostatics derivation of the energy , linear momentum , and intrinsic angular 
momentum equations for an isotropic photon gas . Interactions between the photon 
gas and a medium will be considered only in a general manner. The equations 
obtained from the photon gas formulation are shown to correspond to an electro- 
magnetic formulation of transverse waves. In turn, the direct quantization of 
the electromagnetic variables describing transverse waves are related to the 
phase space representation of the photon gas . The isotropic photon gas relations 
are shown to be consistent with a relativistic Lorentz transformation of the four- 
stress tensor; thus, the relativistic nature of the photon gas temperature, by 
which the photon gas entropy is established as a relativistic invariant , is provided . 

A brief treatment of a nonisotropic photon gas is given as an example of 
the results produced by this nonequilibrium statistical thermodynamics approach. 
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Studies of radiation have contributed substantially to the development of basic 
physics in electromagnetics , thermostatics , statistical mechanics , quantum mechan- 
ics , and relativity . It now appears that the study of radiation also may aid in the 
development of nonequilibrium thermodynamics . 

The generous assistance by and gracious support of the National Aerospace 
Laboratory of Tokyo, Japan, during a significant portion of work on the report 
while the author was in Tokyo is hereby acknowledged . In particular , the author 
would like to thank A. Tate of the National Aerospace Laboratory for his numerical 
computation of the ellipsoidal integrals. 


SYMBOLS 


A 

A .. A .. 
ni] , nji 


A 

V 



electromagnetic vector potential (eq. (Cl)) 

emission, absorption cross section for species n, between energy 
states i and j 

arbitrary four-dimensional vectors (appendix G) 


a 

a ,a 

y V y V 


B 

n 

B 

B 

y V 



c 


cn 

cfj 

n 

D 

D 

E 

E 


arbitrary constant (appendix F) 

simple Lorentz transformation matrices (eqs. (G7) and (G8)) 

scattering cross section for species n 

magnetic induction vector (eq. (B2)) 
arbitrary four-dimensional tensors (appendix G) 

speed of light 
frequency 

photon velocity vector 

electric displacement vector (eq. (Bl)) 
interaction operator (eq. (4)) 
electric field 
electric field vector 
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A 

e 


F 

f 

fP 


g 


ni 


g 

H 

H 

h,h 

hen 

hn 

I 

n 

k 


k ,k® 
n n 


L 

A 

£ 

M 

A A 

m 


m 


N 


n 


n 


P 


A 

P 


energy density (table I), exponential constant (eq. (30)) 
unit vector in polarization direction (eq. (C2)) 
defined mathematical function (eq. (Ell)) 
photon distribution function in phase space (eq. (2)) 

Planck photon distribution function 
degeneracy of ith energy level of species n 

photon gas momentum (table 1) , electromagnetic momentum (eq. (B14)) 

magnetic field 

magnetic field vector 

Planck constant, h = h/2Tr 

energy 

momentum 

specific intensity (eq. (8)) 

Boltzmann constant (eq. (25)), gray gas absorption coefficient 
(eq. (A15)) 

relative , absolute absorption coefficients 

number of dependent variables 

intrinsic angular momentum density (table I) 

magnetization vector <eq . (B2) ) 

intrinsic angular momentum density flux (table 1) 

intrinsic angular momentum number density (eq. (14)) 

molecular particle number density for species n 

photon number density (table I) , dummy numeral (appendix F) 
total-pressure tensor (table I) 
polarization vector (eq. (Bl)) 
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p 

A A 

P 

q 

q- 

r 

S 

S 

s 

A A 

T 

T 

t 

U 

u 


A 

V 

V 

y 

V 


V 


AAA 

W 

X 

y 

X 


x.y ,z 


a 

a 


A A 

3 

Y 


Static pressure (eq . (44) ) 
static-pressure tensor (eq. (H32)) 
energy density flux vector (table I) 
coefficient of infinite series (eq. (Cl)) 

position vector in physical space (x,y,z) 

Poynting vector (table II) 

total entropy (eq. (47)) 

entropy density , length in four space 

Maxwell stress tensor (table II), tensor temperature (eq. (HIO)) 

temperature 

time 

diffusion velocity (eq. (17)) 
dummy variable (appendix F) 
average velocity (table I, eq. (H29)) 
four-space velocity (eq. (G3)) 

volume 

charged-particle average velocity (eq. (B6)) 
stress flux triad (eq. (H24)) 
four-space position vector (eq. (Gl)) 

dummy variable (appendix F) 
physical space Cartesian coordinates 
acceleration vector (eq. (6)) 
defined quantity (eq. (H6)) 

Lagrangian multiplier tensor (eq. (H4)) 
relativistic contraction factor (appendix E) 
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S3> -0 -e- o Q> q A> o> 


unit tensor 


6 

6 relative relativistic contraction factor (eq. (55)) , charge density 

(eq. (B 3) ), variational operator (eq. (28)) 

6 , Dirac delta function 

Ap 

, n wave-number vector , magnitude 

polar-angle coordinate (eq. (H8)) 

,K.. wave-number vector , magnitude (eq. (C2)) 

A A 

wavelength 

. Lagrangian multiplier (eq. (28)) 

distance along intensity direction (eq. (9)) 

Stefan-Boltzmann constant (eq. (41)) 
entropy density flux (table I) 

^ stress-energy density four tensor (eq. (KIO)) 

elliptic double integral function (eq. (H16) , table H-1) 
arbitrary function (eqs. (11) and (12), table I) 

spherical-angle coordinate (eq. (H8)) , dummy variable (appendix F) 
unit direction vector 


Superscripts: 

0 rest condition 

+ ,- positive, negative angular momentum 

' alternate , moving system 

* complex conjugate 
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Subscripts; 

i.j.k 

m 

n 


T 

x,y ,z 


n 

X 


y ,v,(p,x 
0 




dummy indices , energy state i ,j ,k 

medium 

species n 

total , rest 

Cartesian components of vector 
wave-number space variable or operator 
per wavelength 

dummy indices for four-dimensional tensors (appendix G) 
rest condition , transverse electromagnetic wave 
positive , negative angular momentum 


Mathematical notations: 

3 

d r physical space volume element 

3 

d n wave-number space volume element 

" vector , first-order tensor 

"tensor," second-order tensor 
""" triad, third-order tensor 

< > wave-number space average (eq. (11)), macroscopic value of 

electromagnetic variable 

GENERAL FORMULATION 


In this section , the Boltzmann equation for photons , the radiation transfer 
equation , and the equations of change are presented . 


Boltzmann Equation for Photons 

The quantization of a transverse electromagnetic field is obtained by 
direct analogy with the quantum mechanics of particle motion (refs. 14 to 16) . 
Maxwell’s equations for transverse waves are written in Hamiltonian form with 



generalized momenta and coordinates . From this form , the analogy with particle 
motion leads to a corresponding Schrddinger equation for the transverse field. 

The solutions of this equation can be put in the form of characteristic energy or 
momentum functions with photons as a physical interpretation . 

To begin a photon gas formulation , it is sufficient to accept the validity of 
Planck's hypothesis that the momentum of a photon may be considered to be propor- 
tional to a vector wave number fj . The vector velocity of a photon cfi/n . where 
c is the speed of light and n is the vector magnitude , leads to the relativistic 
relation momentum times velocity equals energy . 


hn *c— = hen 
n 


( 1 ) 


where h is the Planck constant . The statement that a photon has a physical posi- 
tion denoted by r and a wave number denoted by n will not , by definition , imply 
a specification that is greater than that allowed by Heisenberg's uncertainty princi- 
ple. Similarly, the statement that a photon has a frequency cn at a time t will 
imply a specification that is less than the maximum allowed by the uncertainty prin- 
ciple . With these implications , the translational state of a photon at a time t is 
represented by a wave number fi at the position r. A specification of the spin 
orientation is also necessary since one of two independent spin states exists for 
each photon . Because the spin results in a photon having an angular momentum 

A 

(ref. 14) of ±h^, the spin state is denoted by + or - , respectively. 

For convenience , the translational state of a photon is represented at a time 
t as a location in a six-dimensional phase space, or y -space, having three physi- 
cal coordinates and three wave-number coordinates . For simplicity , orthogonal 
Cartesian coordinates r = (x,y,z) are used for physical space, and a locally par- 
allel wave-number space fi = (n ,n ,n ) is defined. The motion of a photon in 

X y z 

time is represented by a path or trajectory in this y -space; and, at a given time, 
each photon in the radiation field is represented by a point in this y -space. Be- 
cause photons are indistinguishable , this designation does not represent a unique 
specification of the radiation field; any number of photon locations can be inter- 
changed without a change in the y -space representation of the photon gas . 

The distribution of the photons in the defined y -space is represented by 

densities or distribution functions f^(f,fi,t) and f (f,fi,t). These distribution 
functions integrated over y -space give the total number of photons that make up 
the radiation field. More explicitly, the number of photons in the y -space volume 

element d'^y d'^r about the point (r,fi) at time t is the following. 


[f^(r,n,t) + f (r,n,t)]d^n d^r (2) 
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The total number of photons in the gas at this time is given by 


y' d^r[f’*’(f,n,t) +f (r.n.t)] (3) 

fi-space r-space 


The Boltzmann equation in phase space is analogous to a continuity equation 
in physical space . The rates of change of the photon distribution functions are 
considered for a differential element in y -space and in time . However , the 
resulting analysis is valid only when the changes are considered for a finite volume 
in y -space and for a finite time increment . 

The distribution functions are explicit functions of position , wave number , 
and time . In general , the behavior of a photon gas will also implicitly depend on 
parameters that describe the interactions of photons with material particles. It 
is assumed that changes in the distribution function arising from the statistical 
effects of interactions are equal to those arising from explicit changes in position , 
wave number, and time. Because the Boltzmann equation is written as a rate of 
change with time , this assumption neglects the interaction time with respect to 
the time between collisions (ref. 17) . 

The Boltzmann equations for the rate of change of the photon distribution 
functions are given by 



dr 

" ar" 


dfj 

dt 


v^)f*(f,n,t) = D*f-(f,fi,t) 


(4) 


where D is the interaction operator . The two equations (positive and negative 
polarization) in equation (4) are written for a fixed volume element in y -space 
during the time dt . The first term on the left side is an explicit rate of change , 
whereas the second and third terms are the changes beeause of a ehanging flux 
in configuration space and wave-number space, respectively. The right side of 
equation (4) is an operator representation of the rate of change (due to interactions, 
i.e. , absorption, emission, and scattering) of the distribution function at f, 
n , and t. More detailed consideration of this interaction term is provided by St. 
Pierre (ref. 18) and by others (refs. 14 to 16, 19, and 20) . An illustration of 
the interaction operation between a photon gas and a physical gas is given in ap- 
pendix A . If the two Boltzmann equations contained in equation (4) are added , 

a single Boltzmann equation operating on f^(f,fi,t) + f (r,n,t) is obtained. 

The velocity of a photon in configuration space cfi/ri is equal to the rate 
of change of photon position with time . That is , 


dr _ 
dt = 



n 


(5) 
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The velocity of a photon in wave -number space is proportional to the photon accel- 
eration in configuration space and is defined as 


Id - 

dt ~ 


o_ 

c 


A 

X n 


( 6 ) 


This definition is consistent with relativistic requirements (refs. 21 and 22) and 

retains acceleration vector a“(r,t) with the dimensions of an acceleration or force 

per unit mass . This definition of a" is meant merely to be representative and 
not inclusive of the types of external forces that might be included in the specifica- 
tion of photon behavior . In the absence of an external force or its equivalent , 

and a are set equal to zero. 

By using the expressions for the photon velocity and accelerations , the 
Boltzmann equations are written as 


3 t 


+ 




(f,f|,t) = IiV(r,fi,t) 


(7) 


If the effects of interactions on the photon distribution functions are negligible , 
then the right-hand sides (positive and negative) of equation (7) are set equal 
to zero. For this case, the distribution functions depend only on the boundary 
conditions of the y -space and on the initial conditions. 


Radiation Transfer Equation 

The radiation transfer equation can be obtained directly from the Boltzmann 
equation merely by a change in definition (ref. 2) . Instead of the photon distri- 
bution functions, the intensity of radiation at a given frequency (or wave-number 
magnitude) traveling in a specific direction is considered. The specific intensity 
I is defined as a scalar energy flux of a given wave number , per steradian , at 

A 

r , at time t , and in a direction denoted by U . Thus , the wave-number vector 
is replaced by its magnitude and by the unit direction vector 5 , and the specific 
intensity is written as 

I^(r,n.fiit) = hc^n[f^(r,Ti ,n,t) + f (r,n.i^.t)J (8) 
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As is customarily done , the two orthogonal spin states are neglected in this intensity 
definition . The physical distance along the direction n is denoted by 5 . Then , 
the second operation appearing in the Boltzmann equation is transformed as follows . 



(9) 


For illustration , a and a are assumed to be zero and n is neither a 
function of 5 nor an explicit function of time . Thus , adding the Boltzmann equa- 
tions (eq. (4)) and multiplying by hen yields 


1 J_ 
c at 


^ (? . n ,t) 


( 10 ) 


The prime on the interaction operator denotes the operation on the specific intensity 
as opposed to an operation on the distribution functions . 

The important philosophical aspect of equation (10) is that it represents 
a division of a vector energy flux into an essentially infinite number of scalar spe- 
cific intensities in an infinite number of directions . The equation has a general 
solution for the specific intensity (ref. 1); however, this is not a specification 
of the radiation field as a whole. The practical application of equation (10) (and 
its solution) has been limited to one-dimensional problems . To obtain a solution 
that can satisfy general boundary conditions is difficult , primarily because only 
one differential equation exists in seven independent variables: three position 
variables , three wave-number variables , and time . A useful analogy for develop- 
ing potentially more tractable methods of handling radiation problems is obtained 
with the techniques used in molecular gas dynamics — hence , the concept of the 
photon gas . 


Equations of Change 

In this section , the relation between y -space variables and macroscopic 
parameters in physical space is outlined. The general equations of change and 
the special cases of energy, linear momentum, and angular momentum conservation 

are derived from the Boltzmann equations . The incompleteness of these equations 
in the general case is also discussed (ref. 19) . 

The study of a photon gas as an entity involves consideration of average 
quantities such as photon gas energy , linear momentum , and angular momentum . 
These quantities are obtained by integration over wave-number space to obtain 
average values in physical space . The average values contain an inherent weighting 
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in wave-number space based on the actual distribution functions . Thus , the 
y -space quantity (j>^ has an average value defined by 




The average value for the whole photon gas is given by 

<(|)(f,t)> = <(j)^(f,t)> + <(}> fr,t)> 


( 11 ) 


( 12 ) 


Physical space variables or average values of interest are defined here as 
given in table I . Appendix B contains the corresponding electromagnetic trans- 
verse wave physical space densities , which are compared to the photon gas quan- 
tities in table II (refs. 19, 20, 23, and 24) . A verification of this comparison is 
given in appendix C , in which the electromagnetic quantities are quantized to pro- 
duce the corresponding photon gas relations for energy momentum and pressure 
(ref. 14) . 

It should be noted that 


trace 



+ 



(13) 


A, Ar 

where P is the total-pressure tensor and e.^ is the total energy density, that 


trace (m) = m = hc(n^ - n) 


(14) 


where m is the intrinsic angular momentum density flux, m is the intrinsic angu- 
lar momentum number density , and n is the photon number density , and that 


i = -Cn'^V'^ - n~V~) 


(15) 


A 

where Jl is the intrinsic angular momentum density, h = h/2Tr, and V is average 
velocity . Defining net convective velocity 
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V E |(v'" + v) 


(16) 



and diffusion velocity 


then 


u = 


- 


(17) 


n V 


+ n V = nV 


+ 

he 


(18) 


£ = ilv + 


(19) 


To obtain equations of change, the Boltzmann equations (eq. (7)) are oper- 
ated on from the left with cj)7 and integrated over the entire wave-number space. 

By employing partial integration , by considering c, n , and t to be independent , 
and by assuming that any physical quantity will have a regular behavior at the 
limits of an infinite wave number, the equations of change become 


3 + r it 

-T-T- < <p > + V • <c — d) > 

n n 


( 20 ) 


f+ +.+ 3 



n + 

C—‘V 6 

n 


a 


• V d) ) 

n n 


The general equation of change is the sum of the positive and negative cases of 
equation (20) . 


_3_ 

3t 


< (|) > 


+ 


• <c^- 

n 


n 




+ 



3t 


+ 


C^. 



+ 



' V 

n n 


+ 


3<i 

1 

3t 


+ c3-*Vd) + 

n n 



X 


• V d) > 

n n 


( 21 ) 
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TABLE I.- PHYSICAL SPACE DENSITIES FOR THE PHOTON GAS 



+ 

Photon gas <|(> > 


Number density 

n 

+ 

n 

1 

Number flux 

n V + n V 


c^ 

n 

Total energy 
density 

®T 

+ 

®T 

hen 

Momentum density 

S 

+ 

g 

A 

hn 

Energy density flux 

q 

" + 
q 

V, 2- 

he n 

Intrinsic angular 
momentum density 

s. 

"+ 

±71^ 

n 

Total-pressure tensor 

p* 

P“ 

hcoi 

n 

Intrinsic angular 

momentum density flux 

m 

A ^ + 

m 

A A 

.TlcHi 

n 

Entropy density 

s 

+ 

s 


Entropy density 
flux 

a 

+ 

a 







TABLE II.- COMPARISON OF PHOTON GAS AND ELECTROMAGNETIC 


PHYSICAL SPACE DENSITIES® 


Photon gas 


Total energy density: 


-JhcT]{f'^ + f )d^n 


Energy flux density: 


= + f )d^n 


Momentum: 


-f. 


g -lhr]{f^ + f )d^n 


Transverse electromagnetic wave 
(Gaussian units)^ 


Total energy density: 




Poynting vector: 


<S> = X H> 

4t7 


Momentum: 


<go> 


<S> ^ 1 

2 4ttc 


<E X H> 


Pressure tensor: 


V =/hc!^(f'^ + f )d^n 

J n 


trace ( P ) = e„ 


Pressure = minus Maxwell stress tensor: 


^ ^ 1 9 9 ^ ^ 

-<T > - 5^<E + H >6 

0 8tt 


^<EE + HH> 
4it 


trace 




<eo> 


Reference 19. 


E = electric field vector; H = magnetic field vector. 
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The energy , linear momentum , and intrinsic angular momentum equations for a 
photon gas are obtained immediately from equation (21) , 

Energy: 



V.q = 



( 22 ) 


Linear momentum: 


-|f + v-P = Jhn(3Y + D f)d^ 
n 

+ ^ S~) C23) 


Angular momentum: 


ii 

at 


+ 


V *m = 



+ 


1 (^+ 
— \a 
c ^ 


X 



(24) 


In equation (22) , q is the energy density flux vector; in equation (23) , g is 
the momentum density . Equation (24) is the total angular momentum equation 
minus f cross the linear momentum equation since the linear momentum equation 
is satisfied independently . The photon gas with intrinsic angular momenta is an 
explicit example of a Cosserat continuum (ref. 25) . 

Equations (22) to (24) are the conservation equations for a photon gas. 
Appendix D contains a comparison of the equations of change obtained from either 
the photon gas formulation or the electromagnetic derivation of appendix B . (See 
also refs. 19 and 20.) The general equation of change, equation (21) , can be 
used to generate higher order moment equations or constitutive equations such 
as a momentum flux or stress equation or an angular momentum flux or torsion 
equation. These higher order moment equations will not be required in the cur- 
rent formulation. 

Equations generated from the equation of change , such as conservation equa- 
tions (22) to (24) , are never complete. There are always more dependent variables 
than equations . For example , the 7 scalar conservation equations are written here 
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in terms of 28 scalar components of e„ , P^, q, g, i, and m. The 10 existing 

algebraic relations (i.e. , q = c g, = trace (P) , P, and m are symmetric) 

reduce the number of unknown variables to 18. To form a complete set of conserva- 
tion or momentum equations requires additional information , independent of the 
information obtained from the general equation of change . In this formulation , 
additional information is obtained from the variational principle of statistical thermo- 
dynamics as provided in the next section . 


STATISTICAL THERMODYNAMICS 


In this section , the most probable distribution function , the physical identifi- 
cation of independent variables , the complete conservation equations , and limiting 
cases are presented. 


Most Probable Distribution Function 

The basic principle of statistical thermodynamics is the assumption that 
the entropy of a system is a maximum subject to the independent constraints on 
the system (refs. 7 to 13 and 26) . Each independent constraint is statistically 
weighted by a Lagrangian parameter, which then becomes a dependent variable 
in the equations of change. The most probable distribution function, producing 
the highest entropy, generally is a function of all these dependent variables. 

In principle, the "average values" in physical space can be expressed in terms 
of these dependent variables; therefore, the number of independent equations of 
change is equal to the number of dependent variables . By this procedure , a com- 
plete set of equations for the system is produced . If the state of a system at a given 
point in space and time can be represented by a specified number of dependent 
variables at that same point, then this thermodynamic principle can be applied to a 
single-point distribution function as is carried out here . On the other hand , as 
pointed out by Kogan (ref. 27), in a rarefied situation, the characteristics of a sys- 
tem at a point are highly dependent on the characteristics of the system at neigh- 
boring points . To apply statistical thermodynamics to a rarefied problem requires 
either consideration of a higher order , multipoint distribution function or acceptance 
of the dependent variables in a single-point distribution function as being repre- 
sentative averages over a physical space volume or a time period (or both) . There- 
fore , the results of this formulation are not directly applicable to highly rarefied 
or, in the case of radiation, optically thin situations. The degree of rarefaction 
for photons must be measured by the interaction operator D , since photon-photon 
interaction cross sections generally are negligible. 

For illustration , photons having positive angular momentum +h3 will be 
considered first and the more general case later . Since photons obey Bose-Einstein 
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statistics , the entropy density is , as given in table I , 


+ 


s 





( 25 ) 


where k is the Boltzmann constant . The entropy will be largest if the distribution 
function is the most probable distribution function possible for the constraints 
on the system. In this formulation, all constraints are of the form of "average 
values" at a point in physical space and time. 



(26) 


where i indicates one of a series of constraints . In general , the constraints could 
be statistical averages over a smaller or greater portion of the y -space. For exam- 
ple , in classical thermostatics , the constraints on the entropy are of the form 



'i 

r n 


(27) 


If the constraints are restricted to the form of equation (26) , the condition 
that the entropy density is a maximum subject to a given series of i constraints 
can be expressed as the variational equation 



= 0 


(28) 


where the variables h are the parameters weighting the constraints and serve 

^ + 

as Lagrangian multipliers . If the functions (j) are not explicit functions 
+ + + 

of f , Vf , forth, the Euler-Lagrange equation (eq. (28)) is simply 



In (l 



+ f^ In 



i%i 


= 0 


(29) 
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This equation is readily solved for the most probable distribution function , 




+ + 






- 1 


( 30 ) 


This most probable distribution function can be substituted back into average val- 
ue equation (26) to obtain c|)^ [Xj (r,t)] so that the functions X. become the 
dependent variables for the equation of change. 


With a given number , say L , of independent constraints of the form of equa- 
tion (26), there are, in turn, L Lagrangian multipliers or dependent variables. 
The entropy of the photon gas is a function of all the constraints on the system or, 

in general, all the L dependent variables s^(X^) . The total differential of the 
entropy may be written 


ds"^ = 




(31) 


Incorporating equations (25) and (30) , 


ds = kX 



(32) 


Recognizing that equation (32) contains the differential of equation (26) , 


ds^ = kXj d (j>^> 


(33) 


This , then , is the differential expression between thermodynamic variables that 
provides the thermodynamic significance of the Lagrangian multipliers for each 
value of i . 



(34) 
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Physical Identification of Dependent Variables 

The physical identification of the Xj requires a selection of particular phys- 
ical constraints. In this formulation, the constraints are energy and momentum. 


d)^ = he n 

n,- 



The distribution function (eq. (30)) becomes 


4 * + 

x,hcn+ x„ 


hh 



(35) 


(36) 


Because of the form of the distribution function , it may be readily shown that 


= -n 


Ic—i d n 


^ + 


(37) 


The distribution function becomes 


f^ = 


+, [. fi 

XjhcnU— •- 


-1 


(38) 


If is zero , then X^ is readily identified as 1/kT^ when T^ is the tempera- 
ture describing the familiar Planck distribution function. Since these parameters 
are independent 


X 


+ 

1 



(39) 
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The distribution function for the case of energy and linear momentum con- 
straints can be written as 


f 


+ 



( 40 ) 


Incorporating the Stefan-Boltzmann constant , 


a 


O 5. 4 
2tt k 

icv,3 2 

15h c 


( 41 ) 


the average physical values of interest can be obtained in terms of temperature 
and velocities as given in appendix E . Integrals required for obtaining the relations 
given in table II are given in appendix F . 

The total energy densities can be written as the sum of static energy densities 
and kinetic energy densities 


+ 



+ /^ + ^ + 


e + g *V 


where 


e 





( 42 ) 


( 43 ) 


is the static energy density and y is a relativistic contraction factor. Also, the 
total-pressure tensor P is composed of two parts , a static pressure p and a 
momentum flux 


P = p 6 + g V 


( 44 ) 


where 6 is a unit tensor and 


D- = 2o 4„ 4 _ e 

P - 3c^± - ~T 


( 45 ) 
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The relationships thus obtained are consistent with a simple Lorentz transforma- 
tion (appendix G and refs . 21 and 28) , provided that the static temperature is 
considered to be a relativistically contracted total temperature (ref. 29); that is, 


(46) 

"T “ ~ 


Equation (46) is consistent with the invariant or true scalar characteristics of 
the total entropy 





(47) 


where the rest or total volumes are given by 



= w V 
' + + 


(48) 


Thus, equation (47) , being only a function of the rest conditions and independent 
of relative motion, is a true constant. 

The average values of energy and momentum for a point in physical space 
and time are sufficient to complete the conservation Equations for a photon gas . 
Appendix H contains an example of a more general constraint in which momentum 
flux or stress is the quantity of interest as opposed to energy. The relationships 
obtained for a nonisotropic photon gas are obvious generalizations of the isotropic 
case through extension of the scalar temperature concept to one of a second-order 
tensor . Conceptually , this generalization is immediately obvious . The scalar tem- 
perature is a statistical weighting parameter for the "conserved" quantity energy . 
A complex system of a large number of particles has many more degrees of freedom 
than are accounted for by the single-particle degrees of energy, linear momentum, 
angular momentum , and mass . The anisotropy in a photon gas can be produced 
by an anisotropy in the physical media or by constraining boundary conditions . 

In general , this anisotropy is distinct from relativistic effects alone . 


Complete Conservation Equations 

The two independent spins of a photon give rise to two independent photon 
gases . Each gas has an energy equation and a momentum equation . These equations 
can be coupled to the extent that the physical medium creates or destroys photons 
of a given spin. It is not possible, however, to obtain a single distribution func- 
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tion that is descriptive of a photon gas having a given amount of energy , linear 
momentum , and angular momentum . Both distribution functions are required for 
this description. A photon gas with energy, linear momentum, and angular momen- 
tum can be easily obtained from a combination of the preceding statistical thermo- 
dynamics results . The optimization relations for the + and - photon gases are 
readily combined if the two relative volumes are coincident . 


= V_ = V 


(49) 


In addition , if only one energy constraint is imposed , then the multipliers for energy 
must be identical . 


= T_ = T 


(50) 


The average quantities of interest for a net photon gas are readily obtained and 
are shown in appendix I . 

The familiar thermodynamic relationships hold for the component and for 
the net photon gas . As an example , for the net gas , 

e = sT - p (51) 


An alternative form of the average quantities of interest is obtained by incorporating 
a convective velocity V and a relative diffusion velocity U as in equations (16) 
and (17). 


V = + vj (52) 

U 5 - V_') (53) 
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For convenience, a contraction and^a relative contraction, respectively, are defined 
in terms of the velocities V and U . 


2 


y 


(‘ - 


« 2 2 
2y+ Y_ 



6 



2V*U 2 
- 2 -^ 


c 



(54) 


(55) 


The relationships in appendix I can be rewritten in terms of V and U , as shown 
in appendix J. 


The total-pressure tensor P contains a hydrostatic component , a linear 
momentum flux , and an angular momentum diffusion flux . 


P = + g(v + 6^u) + 


(56) 


The trace of equation (56) gives the total energy density e^, which also has three 
parts . 


e,j, = 3p + g-(v + 6^u) + U-2 


(57) 
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The term -g- U is equivalent to an angular velocity . 

The seven scalar equations of motion, equations (22) to (24) , can now be 
written in terms of the seven scalar dependent variables by the use of appendix E . 
The external forces and the interaction operation (e.g. , appendix A) also must 
be specified. In principle, the solution of a problem also requires coupling with 
the energy , linear momentum , and angular momentum equations of the medium in 
which the photon gas exists . 

A complete set of equations of change for an isotropic photon gas now has 
been established. As shown in appendix G, these equations of change are con- 
sistent with relativistic requirements . To show this consistency for the angular 
momentum requires consideration of £ as the antisymmetric part of a second-order 
tensor (ref. 21). The equations of change are too complex to allow analytical 
solution other than in simple special cases . The following section on limiting cases 
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illustrates some of the simplifying assumptions that make the equations of change 
for an isotropic photon gas considerably more tractable . 


Limiting Cases 

Although laser radiation and many astrophysical sources of emission can 
be highly polarized , thermal radiation is often close to a random polarization state . 

'' 2 

If a random polarization is assumed, then U and 6 are zero, as is the solution 
of the angular momentum equation (eq. (24)). For this case, the physical space 
dependent variables take on a much simpler form , as given in appendix K . 

The limit of V/c approaching unity corresponds to a unidirectional photon 
flux in which temperature, pressure, and entropy approach zero. In this limit. 


lim g 




2 


c 


(58) 


and, for a finite momentum, 


lim 


V 

c 






constant 


(59) 


At the other extreme, if V/c is small compared to 1, the binomial expan- 
sion of the variables in appendix K is useful as given in appendix L . This approxi- 
mation provides for less cumbersome variables and should be applicable in media 
for which temperature gradients are not extremely severe. 

Although the equations are complete in principle , analytical solutions are 
difficult except for the most elementary cases . In addition , the practical application 
of the photon gas formulation to situations involving finite boundaries or walls 
requires a specification of physically realistic boundary conditions. At a wall, 
however , an optically thick restriction generally is not valid . Thus , the applica- 
tion of this photon gas formulation to finite boundaries requires, in general, 
restricted equations at the wall or , in analogy with physical gas dynamics , a 
boundary-layer or slip-layer formulation. 
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CONCLUSIONS 


A statistical treatment of photons , analogous to molecular gas dynamics , 
has been used to provide equations of motion for a thermal radiation continuum . 

The photon gas equations have been shown to be identical to equations of motion 
obtained from an electromagnetic formulation based on Maxwell's equations written 
in a microscopic form . The correspondence between the two formulations has been 
established by a quantization of the transverse electromagnetic wave form of energy, 
momentum , and stress . Nonequilibrium statistical thermodynamic principles have 
been applied to the photon gas to determine the most probable distribution function 
or the highest possible entropy subject to the constraints on the system. This 
application has provided a truncation of the equations of change or the moment 
equations at an arbitrary level. The maximum entropy principle has been applied 
to an equilibrium statistical thermostatics derivation of the energy, linear momentum, 
and intrinsic angular momentum equations for an isotropic photon gas. These 
conservation equations are obtained in terms of a temperature , an energy average 
velocity, and an intrinsic angular momentum diffusion velocity. Comparison of 
these relationships with a Lorentz transformation has provided the relativistic 
variation of the isotropic photon gas temperature and established the entropy as 
a true scalar. 

An example of extending statisticeil thermodynamics to a more general non- 
equilibrium application such as a stress or momentum flux constraint has also been 
presented. This process results in the definition of a second-order tensor tempera- 
ture that is descriptive of the nonisotropic gas or of a stress state. 

Although this photon gas formulation provides , in principle , a complete 
set of differential equations , the solution of these equations for a practical problem 
is quite complex. In addition, the treatment of more general nonequilibrium problems 
such as radiation in optically thin media, the boundary equations at a medium inter- 
face , or the spectral distribution of radiation in a real medium requires extensions 
of this photon gas formulation by consideration of more general constraints . 

The prime application of this work resides in rendering an improved descrip- 
tion of thermal radiation transport; potentially, an important application also exists 
in producing a vivid demonstration of the power of statistical thermodynamics from 
the vantage point of information theory. Information theory permits a rational 
truncation of statistical moment equations based on probability concepts rather 
than on simple particle mechanics . Thus , statistical thermodynamics is applicable 
to nonequilibrium phenomena and to macroscopic phenomena such as fluid 
turbulence . 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, February 18, 1975 
909-44-39-00-72 
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APPENDIX A 

THE INTERACTION OPERATOR 


The right-hand side of the Boltzmann equation represents the net rate of 

photon production due to interactions with the medium . Thus, Z)“f“(r , f) ,t) is 
representative of this net production of photons at the phase space point r , fi in 
3 3 

the volume d r d n and between the times of t and t + dt. The net production 
is the difference between the total production and the depletion or extinction . 

The total production consists of spontaneous emission, induced emission, and scat- 
tering into the six-dimensional phase space, or y -space, volume element. The 
depletion is due to absorption and scattering from the y -space volume element. 

The relative importance of various interaction processes depends on the 
physical situation. A comprehensive treatment of photon-media interactions is 
beyond the scope of appendix A; however, the basic aspects of absorption, emis- 
sion, and scattering are considered to the extent that these interactions tend to 
drive the photon gas to a thermodynamic equilibrium with the media (ref. 18) . 

Consider a molecular gas of several species n, each having a particle con- 
centration of [N^] . Spontaneous emission by these particles may proceed from 

an energy state i to an energy state j . 


N 


ni 


+ photon (hn)^ 


(Al) 


The discrete energy levels for the particle are obtained from the quantum mechanical 
description of the possible states for the isolated particle. These states may be 
degenerate in terms of the same energy. This degeneracy for the energy state 
i will be denoted by g . . 

.7 &m 

Let the probable spontaneous emission rate per particle of species n , in 
the i energy state, be represented by cA^jj(ri,t) for photons of energy hen. 

The quantity A^^^ has the dimensions of a cross section and is directly propor- 
tional to the Einstein transition probability of spontaneous emission. The cross- 
section dimensions are retained since only the emission in one specific direction 
given by fj is considered. The total probable rate for spontaneous emission at 
f , n , and t is then 


EEEn j(f,t)cA (n,t) 
n 1 j ■’ 


(A2) 
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The summations are only over possible combinations that can produce a photon 
of wave number fi . In addition to this spontaneous emission , there is an induced 
emission that is proportional to the number of photons present . Thus , the total 
probable rate of emission is 


[l + f‘'(?.n.t)]X;EZNni(r.t)cA^ij ffi.t) 
n i j 


(A3) 


The probable rate of absorption per particle of n in the state j will be 
represented by the cross section A^^^(fi,t). This cross section is for the reverse 

of the reaction shown in equation (Al) . The total probable absorption rate at 
r , fi , and t is 




n 1 1 


nji 


(A4) 


The difference of equation (A3) minus equation (A4) is the net rate of photon pro- 
duction due to absorption and emission. 


£>~f" 


(l + f-)N ,A„,. - f-N .A .. 


n i j 

L ' ni nij nj nji 


absorption and 
emission 


(A5) 


If the cross section for a particle of species n to scatter a photon of wave 
number fj to a photon of wave number rl' is defined as B (fi,fi'), then the net 
rate of photon production due to scattering is 


D-f- 


n 


d n 


B (f)' 

n' ' 


,t)N^(r,t) - B^(fi 




n.t) Nj^(r,t) 


(A6) 


scattering 


If the molecular gas is assumed to be in a state of local thermodynamic equi- 
librium at a local temperature T(r,t) , then, by Boltzmann statistics, 




^ni 


g, 


n] 


hen 

kT 

e 


(A7) 


28 


where k is the Boltzmann constant . The energy hcri is the difference between 
states i and j . In addition , if the gas is in equilibrium , the principle of detailed 
balance requires that the probability for a microscopic process must equal the 
probability for the reverse process . This relationship provides the rationale for 
the assumption that 


g .A .. = g -A .. 
^ni nij ^nj nji 


(A8) 


The net rate of photon production for absorption and emission may now be 
written 


D-f- 


n i j 


nj nji 


,(l 


absorption and 
emission 


+ 


hCn 



(A9) 


or 


D-f- 


= ck 


0 


hen 

kT 


absorption and 
emission 



(AlO) 


A true absorption coefficient has been defined as 

(n.t) (All) 

n i j 


The induced emission term appearing in equation (AlO) appears to be more 
like a negative absorption . Thus , an effective absorption coefficient is defined 
as 



(A12) 
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and a Planck photon distribution function as 


f^[T(r,t)] 



(A13) 


These two definitions may be incorporated into equation (A 10) to obtain 


±,.± 

D f 


-ck (f- 


absorption and 
emission 



(A14) 


The physical significance of r may be obtained readily by consideration 
of the case in which the solution of equation (A14) is zero for an arbitrary absorp- 
tion coefficient . This value implies that the molecular gas and the photon gas are 
in a state of mutual or complete thermodynamic equilibrium . Thus , the Planck 
distribution function is the one that would not be altered by the molecular gas . 

It might be considered as being defined by the temperature T , whereas the temper- 
ature T(r,t) is a local property of the molecular gas. 

As a special case , it can be assumed that the effective absorption coefficient 
is independent of fi . 


k^(r,fi,t) — >• k(f,t) (A15) 

In equation (A15) , k is the gray gas absorption coefficient, and this assumption 
is commonly referred to as the gray gas assumption . In general , it is not a very 
realistic assumption; however, it retains the basic interactions of emission and 
absorption and greatly simplifies the equations . 

The physical features of the scattering interaction term are best realized 
by considering the equations of change . The prime effect of the scattering term 
is a dissipation of the net relative momentum between the medium and the photon 
gas. 
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APPENDIX B 

ELECTROMAGNETIC FORMULATION OF RADIATION 


In this appendix, Maxwell's equations, electromagnetic energy, electromag- 
netic momentum , and an electromagnetic phase space are considered . 


MAXWELL'S EQUATIONS 


In comparing the equations for a photon gas with those for electromagnetic 
waves , it must be noted that the two approaches are complementary but distinct 
(ref. 19) . The comparison in this appendix is in terms of macroscopic or statis- 
tical average quantities . Brackets < > are used to denote statistical averages 
of the enclosed quantities. In general, these averages have a variation with position 
and time (r,t) . This average is analogous to the statistical average for photon 
transport as presented in appendix C . 

The macroscopic electric and magnetic field intensities are defined as <E> 
and <H>, respectively. In addition, the displacement vector and the magnetic 
induction vector are defined as <D> and <B>, respectively, in Gaussian units. 


<D> = <E> + 4tt<P> 


(Bl) 


and 


<B> — <H> + 4ti<M> (B2) 

In^ equation (Bl) , <P> is the polarization vector of the medium; in equation (B2) , 
<M> is the medium magnetization vector. In free space , <P> and <M> are zero. 

The familiar form of Maxwell's equations (in Gaussian units) (refs. 23 and 

24) is 


V« < D > = 4tt6 


(B3) 


V<B> = 0 


(B4) 


V X <E> + — <B> 
c 9t 


= 0 


(B5) 
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and 


V X <H> - = ^6v (B6) 


In equations (B3) and (B6) , the current is the charge density 6 times the average 
charge velocity v. 

The macroscopic form of Maxwell's equations cannot be used for a compari- 
son with the photon transport approach. It is necessary to revert to the microscopic 
form of Maxwell's equations, transform this form to equations of motion, and then 
find the rate of change of the average values (ref. 24) . Thus, the macroscopic 
energy and momentum equations for the electromagnetic field cannot , in general , 
be obtained from the macroscopic forms of Maxwell's equations. This restriction 
is probably one of the major sources of difficulty encountered in attempts to include 
absorption and emission coefficients in electromagnetic field approaches to radiation 
transport. The microscopic form of Maxwell's equations is effectively the same 
as the equations obtained when the brackets are removed from all the terms in equa- 
tions (B3) to (B6). 


ELECTROMAGNETIC ENERGY 


The energy density of an electromagnetic field is assumed to have the form 

<e(r,t)> = ~ <E-D + H*B> (B7) 

OTT 

To compare this electromagnetic energy density with the energy density of a photon 
gas , two components 

<e> - <Cq> + <e^^> (B8) 


where 




and 

<e > 5 J-<E*P + H*M> (BIO) 

m 2 
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are defined. In general, the energy density contains more than just a trans- 

verse electromagnetic wave energy . It is assumed that , for a basically neutral 
and dilute gas, <Cq> may be considered as the transverse wave energy of the 

field. On this basis, <Cq> is compared to the total photon energy density e,p. 

If this assumption is not made, it is necessary to divide into two additional 

parts for a comparison . The energy is the field energy due to the polariza- 

tion and magnetization of the material particles. On the average, for a neutral 
gas, <Cj^> is negligible except for localized fluctuations in the immediate vicinity 

of a material particle. In the absence of a medium, is zero, and <Cq> should 

be identically equal to the total photon energy density e,j, (table I) . 

It is assumed that the Poynting vector 

<S> E ^ <E X H> (BID 

4ir 


represents a transverse wave energy flux. Although some question exists as to 
the general interpretation of <S> as a physical flux of energy, the vector <S> 
appears to be a reasonable assumption for a transverse wave field. As with <Cq>, 

it is assumed that the Poynting vector contains only contributions from the trans- 
verse waves , or that the energy flux associated with longitudinal waves is negligi- 
ble. If these assumptions are valid, then the Poynting vector should be compara- 
ble to the photon energy flux vector q . This definition of the Poynting vector 
differs from the conventional use of the "basic" magnetic field B (ref. 23) . 

This difference is made intentionally to show analogy with the photon gas vacuum 
energy flux. 

The energy equation is obtained from a scalar combination of the microscopic 
forms of equations (B5) and (B6) . By employing the microscopic form of equa- 
tions (B3) and (B4) , performing a statistical average, and using the definitions 
(eqs. (B8) and (Bll)) , the following is obtained. 

A. 

^ /s ^ Si P ^ Si ]VT A. 

<6q> + V.<S> = - + H - - <6E-v> (B12) 


Equation (B12) is to be compared to the photon gas energy equation (eq. (22) 
and table II) . In the absence of a medium, the right-hand sides are equal to zero, 
and the quantities on the left-hand sides should have a one-to-one correspondence . 
When a dilute medium is present , this correspondence depends on the assumptions 
noted previously. It is also noted that if photon scattering processes are conserva- 
tive, they should not appear in equation (B12) . However, absorption and emission 
are represented on the right-hand side of equation (B12) . 
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ELECTROMAGNETIC MOMENTUM 


The momentum of an electromagnetic field is assumed to have the form 


<g(r,t)> 


4ttc 


<D X B> 


(B13) 


It is desirable to separate this momentum into two components . 


<g> = —fr <S> + — <E X M + P X 

^2 c 

c 


(B14) 


In the absence of a medium , the momentum of the field is the Poynting vector 
2 

divided by c . 

When the Poynting vector is^examined for a simple plane monochromatic 
transverse wave, it is noted that S is not constant. Rather, the value of S 
oscillates with time and periodically becomes zero. Even this behavior, however, 
is not inconsistent with the similar but separate representation of q for one pho- 
ton . It is recalled that the energ 5 ^ of a photon must also be specified over a finite 
period. The Poynting vector <S> is an average over an assumed large number 
of small transverse waves . Similarly , q is the average energy flux for a large 
number of photons . 

To obtain a momentum equation , the vector product of D with a microscopic 
form of equation (B5) is added to the vector product of B with a microscopic form 
of equation (B6) . Vector identities are used to rearrange the result, which is 
then integrated over wave-number space to obtain 


3 

at 


:g> - 7 


<E (v«e) + H (v.h)> 

X fij + P X (v X e) + M X (v X h)> (B15) 


as 


The Maxwell stress tensor for a pure electromagnetic field has been defined 


<TQ(r,t)> E ^ <EE + HH - J 6 (E + H^)> 


I A A A A AA 

= <EE + HH> - 6 <e«> 
4-rr 0 


(B16) 
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where 6 is the unit second-order tensor . The physical significance of this second- 
order tensor is a stress due to the relative momentum flux of the vacuum compo- 
nents of the field . The transverse portion of the tensor should correspond to the 
negative of the photon gas relative pressure tensor P (table II) . 

In^ further development of the momentum equation, the relationships for 
D and B (eqs. (B3) and (B4)) are included to obtain 


1 

2 

c 



7 * <Tq> = - <6 


(e + I - ft). 

+ <E (vf) + H (7-m)' 


. ,1e « + « ft. 

c 9t c 3t 


(B17) 


Equation (B17) will be taken as the final form of the momentum equation and should 
be compared to the photon gas momentum equation (eq. (23) and table II) . 

As in the energy equations , in the absence of a medium and any external 
forces, the right-hand side of equation (B17) is equal to zero. The right-hand 
side of equation (B17) can be manipulated into many forms. The first term is the 
normal form of the Lorentz force , the second term is the photon scattering , and 
the last term is the emission and absorption . 


AN ELECTROMAGNETIC PHASE SPACE 


The discussion of electromagnetic phenomena was initiated with the macro- 
scopic form of Maxwell's equations. The statistical or bracket average used to 
define the macroscopic values of the fields is a classical sum, or integration, over 
an independent wave-number space . The requirement for this statistical average 
becomes apparent when a comparison is made with the photon gas formulation. 

The general transverse electromagnetic field cannot be discussed without 
considering it as being divided into component waves . This division is physically 
reasonable because the field cannot, in general, be associated with the motion of 
one large charge. Because the component waves can have various directions and 
frequencies , it is convenient to consider a wave-number vector . The analogies 
with the photon gas are straightforward . For example , a wavelength greater than 
the extent of the physical volume element under consideration is difficult to observe . 
However , the wave number , the position , and the time must all be independent . 

It is not certain that E and H are the best parameters to describe the transverse 
waves in this space; however, they do offer advantages and will suffice for this 
discussion . 


The effect of the wave-number space average on the form of Maxwell's equa- 
tions is negligible, since the equations are linear in terms of the fields. However, 
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the energy and momentum equations are not linear in the fields , and their forms are 
definitely affected by the wave-number space average , Explicitly , 


<E> <E> ^ <EE> 


(B18) 


This difference can be illustrated by the following general tensor identity , which 
is derived for an explicit relation between Maxwell's stress tensor and the Poynting 
vector . 






(B19) 


This algebraic tensor equation only holds in what might be termed an electro- 
magnetic phase space. In general, equation (B19) does not hold when written in 
terms of the corresponding statistical or macroscopic variables . In the special case 
in which S is equal to zero , no net energy flux occurs in any direction . For sim- 
plicity , it is assumed that the symmetric tensor T^ has been diagonalized. Under 

these restrictions, equation (B19) produces three relations: 


T 


2 

0 ,xx 




(B20) 


However, because the trace of -T^ is e^, equation (BIO) is identically equal to 

zero , which indicates the absence of radiation . If this physical interpretation of 
the Poynting vector is correct, equation (B20) is not valid for the macroscopic 
quantities . However , the analogies drawn in^this study tend to support the physical 
interpretation given to <S>, and Equation (B19) is also valid for 

the six-dimensional phase space, or jj^-space, quantities associated with the 
photon gas parameters e,^,, q, and P. In general, the electromagnetic phase space 

is a six-dimensional p-space. For example, 

E = E(r,fi',t) (B21) 


The prime on the wave-number vector distinguishes the vector from a direct associ- 
ation with a photon wave number. 

The Poynting vector and other quantities familiar to physical space are not 
necessarily fundamental to the u -space. For example, 

= ~[E{?,n',t) X H(?,n',t)] 
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(B22) 


is a double- valued function of ri' for the two polarization components. 

The most important result of this statistical treatment is that the stress tensor 
in a macroscopic sense cannot, in general, be represented by the microscopic form. 
However , this result does show that the macroscopic electromagnetic formulation 
is incomplete in the same sense in which the photon gas formulation is incomplete 
without the application of statistical thermodynamics . The algebraic relation of 
equation (B19) is valid for the macroscopic quantities for a one-dimensional phys- 
ical space. For this case, equation (B19) is merely a scalar identity. 

Although a classical electromagnetic phase space is somewhat abstract , the 
quantum mechanical concept is straightforward. As shown in appendix C, the 
quantized forms of Oq , "Tq , and S are precisely the same as those of the photon 
gas quantities . 

The comparison of the left-hand sides of the energy equations and the 
momentum equation (table 11) is reassuring, both with respect to the physical 
quantities and with respect to the mathematical operations . A comparison of the 
right-hand sides of these equations requires a more extensive study . 
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APPENDIX C 


QUANTIZATION OF ELECTROMAGNETIC VARIABLES 


Correspondence between the variables and the equations of change for a 
photon gas and a statistical (yet classical) form of the electromagnetic formulation 
is given in appendix B . The statistical average or bracket notation , however , 
is a quantum mechanical concept . To establish the correspondence between the 
bracket average and the photon gas wave-number average requires a quantiza- 
tion of the electromagnetic variables associated with transverse waves as given 
for example by Heitler (ref. 14) , It should be noted that the quantization and 
the statistical averaging are conceptually and mathematically distinct. For a clarifi- 
cation of this point, refer to reference 30. 

The electromagnetic vector potential is expanded by the use of a complex 
series of plane waves . 


A = 




* “ix * "ix^x 


+ q, A. 


(Cl) 


A -f 



2c 




(C2) 


+ 




(C3) 


Here , k, is a wave-number vector for the plane wave , e^ is a unit vector in 

A A 

the direction of the wave polarization, and * designates complex conjugate. 
For transverse waves , 


VA^ = 0 


(C4) 


i . 

which makes the e always perpendicular to 
A 

the normalization over the volume integral. 



Equation (C2) is a result of 


+ A + * 2 

A“»A” dF = 4 irc 6. 
Ay Ay 

F 


(C5) 
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where 6^^ is the Dirac delta function . The scalar coefficients in equation (C 3) 
are normalized to the photon density 


±* ± h ± . 

where the photon number may be equated with the photon gas density 

«-=E 

A 


(C6) 


(C7) 


The electromagnetic field strengths are obtained by the operator representa- 
tion of the vector potential . 


E = 


^ "c If ^ ^ A, + q,A, - q, A 


+ /V 4- +*/N +* 

■ '^X“X ^X “X 


H = V X A 


Ek(«x ^ 


X rx " *x 


* ‘>xk ■' *x) - ‘ix’k ” 


* / . 

^x 


(C8) 


(C9) 


where the scalar coefficient time derivatives are replaced with the operator 
representation 


Tt -‘V'5x Wio) 
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The bracket operation for energy , for momentum , and for the total-pressure ten- 
sor is obtained directly from the operators in equations (C8) and (C9) . 


Energy density: 


<Co> 


= <^(e^ + H^)> 

" s " «r ^ “a") 


= EK'IaAa’Aa + 'Jx'lx 


Momentum density: 


= E h='<x(''x * =/>'on(t* * f)a\ 




E X H 


a.^ a. 

E, XH, .H, 


S XxjqX* 


“ ("x " ''x 


♦/V— / + 

+ X (k^ X a 


Since 


A i A +♦ 


(C12) 


(C13) 


(C14) 
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then 


c 


(4- X#^4. ♦ 

'»x'’x*x'*x i^x 


= Z *'^x("x * "xj 


(C15) 


Total-pressure tensor: 


- <Tq> = <e^j> 6 - <EE + HH> 


j ^ A. A A A A A ^ A ^ \ 

= 'V ' - E ^ Hx»x " »x»x) 


^,±v;± 2 ± ±*i±i±* 

'^x^x^'^x Ox'ix ^^x^x 


2 ± ^* ( t 2^^±-^± 

= (4,0 Je^ej 


= hcKjn^e^e* 


(C17) 


In a similar manner, since e and k are orthogonal; 

A A 


^X^X ^ ^®'^x^x('^ ' ®X®X 2 


(C18) 


41 



Therefore 



Thus, the quantization of the transverse electromagnetic field variables substantiates 
the photon gas representation. 
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APPENDIX D 

COMPARISON OF ENERGY AND MOMENTUM EQUATIONS OBTAINED FROM THE 
PHOTON GAS AND AN ELECTROMAGNETIC FORMULATION 


In this appendix, energy and momentum equations obtained from the photon 
gas and an electromagnetic formulation (ref. 19) are stated for comparison. 


ENERGY 
Photon Gas 




n 


Electromagnetic 


3 

3t 


<e„ 


+ V 


< s > 


-<E 



+ H 


3M ^ 

* 3t 


MOMENTUM 


Photon Gas 



Electromagnetic 


_3_ 

3 t 



- V* 





A A 


+ <E(V*P - 5) + H(V*M)> 


( 22 ) 


(Dl) 


(23) 


(D2) 
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APPENDIX E 

ISOTROPIC PHOTON GAS PHYSICAL SPACE DENSITIES 


In this appendix, terms of isotropic photon gas physical space densities 
are defined on the basis of the photon distribution function in phase space f . 


f 


+ 


hen/. V~ 

c 



-1 


1 . 


(40) 


Photon number: 


+ 


n 


20_o_ T. 3 4 
27 kc ± 


(El) 


Total energy: 


e 


~T 


2 - + 
c ± ± 


V 

8 o „ 4 6 

3 c ~2 

c 


(E2) 


Linear momentum: 


g. 


2 


_ 8 £ T 

3c 


4 6- 


(E3) 


Intrinsic angular momentum: 


j ^ 20 h 3^ 4^ 
■ - 27 ^^^± n 


(E4) 


Total pressure: 


AA 



1 1 


V V 

^ 8 a„ 4 6 

3c^± Y± — 2“ 


(E5) 
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Intrinsic angular momentum flux: 


m 


+ 


-27 + 





(E6) 


Entropy: 


s 


- 

± 3 c V 



(E7) 


Entropy flux: 


o = s 


(E8) 


where 





(E9) 


F(A) 



n=2 



2n + I 


(ElO) 


trace 



0 


(Ell) 
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APPENDIX F 

INTEGRALS REQUIRED FOR FORMULATION 


The evaluation of integrals over wave-number space is simplified by appli- 
cation of the following mathematical relationships . In this appendix , x, u, 
and <p are dummy variables , n is a dummy numeral , and a is an arbitrary 
constant . 



(FI) 


oo 



n=l 


OO 



n=l 


4 

7T 

90 


OO 



n=l 


3 

ir 

25,79436 ... 




nl 


0 


IT 



sin <p d<p 
+ a cos (p ) ^ 


+1 


/- 

(1 + 


du 


-1 


(1 + au) 


(l - 


+ 1 


/ 


u du 


-1 


(l + au) ^ 



(F2) 


(F3) 


(F4) 


(F5) 


(F6) 


(F7) 
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du 


2 


(F8) 


+1 


f- 

J a + 


(1 + au) 


-1 




+ 1 


A 


u du 


-2a 


-1 


(1 t au) (j . „2)' 


Also , the Stefan-Boltzmann constant is defined as 


0 = 


2A^ 

ISh^c^ 


(F9) 


(41) 


An example of the application of these integrals to the evaluation of the photon gas 
densities is provided in appendix H for the nonisotropic photon gas . 
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APPENDIX G 


THE LORENTZ TRANSFORMATION 


For consideration of the relativistic aspects of a photon gas , a four- 
dimensional space of position and time is required. The position vector in this 
space (ref. 21) is 


X 



An element of length in this space is an invariant scalar 


(Gl) 


ds = 


4 


2 

c dt 


dr 


(G2) 


The generalized velocity is a four-space vector 



(G3) 


where V = 
dimensional 


g and Y = 



The velocity is a true vector in a four- 


space; therefore, its magnitude is a true scalar. 


V V 
y y 



-1 


(G4) 
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The stress-energy density tensor in this four- dimensional space for the photon gas 
and transverse electromagnetic waves (refs. 14, 22, and 28) is 




4<eQ> 


This four tensor has an invariant magnitude 


X = Qrr, - trace = 
uu T N / 


To transform four vectors or four tensors^from one coordinate system to another 
(prime ' ) moving at a constant velocity V requires the application of the Lorentz 
transformation matrices (ref. 21) . 


iyV 


O A A 

T a. Y VV 

Y + 1 2 

’ c 


My^ 
' c 


The transformation for four vectors and four tensors is obtained by 


A = a A 
y y V V 


A = a A 
V vy y 


f 

B = a B. .. a.. B = a . B a .. 
yv yifi 4>X AV (j)X i}iy yv vX 


L 



If a photon gas is at a complete thermodynamic equilibrium and at rest in 
a given coordinate system , the stress-energy density four tensor is 


0 


I2.T 

3c 0 



0 

+3 


(GIO) 


To an observer moving with a velocity V with respect to this static photon gas, 
the stress-energy density four tensor is 


T 

yv 


y<)) (|) A A V 


i°-T ^ 

3c ^0 


T d 2W 
■6 - 4y -g- 

c 


•d 2 V 

_i4y - 



+3 + 4y 



(Gll) 


This tensor is identical to the stress-energy density four tensor obtained through 
statistical thermodynamics for the photon gas with both an energy and a momentum 
constraint (appendix K) providing 


T 


0 


yT 


(G12) 


Therefore, T is the apparent or static temperature seen by a general observer, 
whereas T^ is the total or rest temperature of the photon gas. On the other hand, 

the distribution function is not affected by relative motion even though the tempera- 
ture and wave number are . The distribution functions in rest and moving frames , 
respectively, are 


f 


0 



f = 


e 



1-1 


- 1 


(G13) 


(G14) 
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If f = fp and Tq = yT, then 


^0 



(G15) 


which is the correct Doppler relationship between the relative wave-number magni- 
tudes for a given photon (ref. 14) . Thus , the characteristics of a photon gas moving 
with a velocity can be obtained either through statistical thermodynamics with a 
momentum constraint or by a Lorentz transformation of a static photon gas . 
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APPENDIX H 

THE NONISOTROPIC PHOTON GAS 


Statistical thermodynamics is classically restricted to the determination of 
the most probable distribution function, producing the highest entropy, subject 
to the constraint of energy . The extension to constraints involving other conserved 
quantities , such as linear momentum and angular momentum , is straightforward 
as given, for example, in this report. However, since the concepts of statistical 
thermodynamics can be obtained from information theory , the principle should 
be applicable to general independent constraints beyond the simple body parameters 
that are conserved (refs. 7 to 13 and 26) . For example, a photon gas may have 
the highest entropy possible under the constraints of the initial conditions , the 
boundary conditions , or the medium and not be isotropic . This case is considered 
as a demonstration of the concept . For simplicity in notation , a randomly polarized 
photon gas is assumed . The extension to the case of intrinsic angular momentum 
is obvious but cumbersome . 

Consider the constraints on the photon gas to be the total-pressure tensor 
and the momentum density 


y '* . . 
he im 
n 


f 

I d n 


(HI) 


g = 2 yinf d'^n 

n 


(H2) 


where 


f = f E f 


(H3) 


The most probable distribution function under these constraints is obtained 
from the variational equation for the maximum entropy (eq. (28)) . 


6 



( 1 + f ) + f In 



+ 


6 : 


P 


2 + 



(H4) 
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Here , 6 and 3 *V are the tensor and vector multipliers for the total-pressure 
tensor and momentum, respectively. Equation (H4) produces the distribution 
function for a nonisotropic , randomly polarized photon gas with momentum . 



(H5) 


The unknown Lagrangian multipliers are, in principle, obtained through 
equations (HI) and (H2) . In this case, however, it is simpler to first obtain the 
entropy in terms of 3 and V. Defining 


\n c/n 


the entropy integral, based on equations (FI) and (F2) , becomes 


s = 2k 


zj(k * 


-nhcan 2 


d n 


If spherical wave-number coordinates are used 


2tt Tt <» 

r ^ —/ -!/■ sin <p d9 


0 0 


Using equations (F5) , (F3) , and (41) , the entropy becomes 

2tt tt 


= ,Jy I ] f ae 

/"'-I 


dtp 


0 0 


2it it 




dtp' 


0 0 


(H6) 


(H7) 


(H8) 


(H9) 
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To describe the nonisotropic state , a tensor temperature is defined . 


T = 



where 


T*3 



6 


(HIO) 


(HID 


iTl =■ 




(H12) 


The scalar temperature with its classical significance is retained. 


T = I trace (t) 


(H13) 


The entropy for a nonisotropic photon gas may be written 


where 


s 


16 a 4 


Tl $ 



(H14) 


(H15) 


and a dimensionless number that is only a function of the rest principle direction 
temperatures has been defined 





2tc 


% 


/./ 


sin 9 d 9 




(H16) 
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In general , $ is an elliptic integral . This integral has been computed 
numerically and is tabulated for the rest condition together with an approximate 
expression in table H-I. It is obvious from this table that $ is a very weak 
function of the deviations from an isotropic temperature condition . 

The entropy for a nonisotropic situation is always lower than the isotropic 
entropy for the same energy , a condition verified by the fact that 


1 /■' 


AA 

j trace v 

rjj > 

T 


(H17) 


The randomly polarized , nonisotropic photon gas variables can be expressed in 
the integral forms 




a 


2kc 


/i[ln ll + tl 

•'n ■- 


+ f In 




7t 


de 



sin cp dq> 

3 

a 


0 


(H20) 


(H21) 
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g E 2 / hnf d^n = 

J c(hc)‘ 


2% n 




W E 2 




2n % 


f d^n 


n^l’' J J ’> “ 


0 0 


The third-order tensor w appears in the third moment or stress equation, which, 
from equation (21) , is 


A, A. 

+ V* w = 2 Df d^n 


yy ^ ^ ^ ^ A 

for D = D =D and a =a =0. Reductions in w are q (in index notation) . 


W... = W... = W... = Q. 

1)J Jl) )]1 1 
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In addition, the trace of equation (H25) is the energy equation. Noting the pattern 
in equations (H18) to (H24) , the following interrelationships may be obtained. 


where 


„ = 40^ ^ 4 

” 27ck ^ ^ * 

(H27) 

II 

CO 

(H28) 

V 

(H29) 

o = sV 

(H30) 

* 1 3 a 

g = 

^ 3V 

(H31) 

P - Vff = - — = n 

* 4k - P 

(H32) 


is the static-pressure tensor . 

VP 

^ . A. 

- VgV 

+ Wg (H33) 


AA A 

^ 3$ 
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PV 
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^ AA 

pv 

+ Vp 
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The expressions containing derivatives with respect to the Lagrangian multi- 
pliers are, for example, 



3s 

3s 

3s 



'^12 


3s _ 

3s 

3s 

3s 

aA 

33 


'^22 

'^23 


3s 

3s 

3s 


3 031 

^^32 

'^33 


(H34) 


Equations (H33) and (H26) may be used to relate the momentum and velocity. 



It is noted that, since 


S: 


3s 

3 $ 


= -3s 


6: -It (in 4>) = 0 


(H35) 


(H36) 


(H37) 


the thermodynamic relationship between entropy and the local nonisotropic pres- 
sure is 


A AA /aA^ a \ A AA AA 

s = |k$:(P - Vg) = ^kg;p 


(H38) 
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Equation (H36) is readily obtained through the indicated operation on equation 
(H20) . In a similar manner , equation (H37) follows from the same operation on 
equation (H16) . 

The momentum and pressure for a nonisotropic photon gas may be written 
in a more convenient form . 


^ 4g 

AA 

4 

S - 3 

3c^ 

T 

Y 'ST* 


4X V + c 


(in ♦) 

- 


P 


4g 

3c 




T 



P 


9 

TU 


nkT 


9 3 1.4 * 

In $ 

19 ae 


(H39) 


(H40) 


(H41) 


Equation (H25) and the momentum equation form complete set of nine equa- 
tions in terms of the nine dependent variables, V and T (i.e. , T is symmetric) . 
The solution of these equations in general would be quite difficult. The first obvious 
simplification is to approximate g and p by the dominant terms in equations (H39) 
and (H40) . 


S 




Q ^ ^ 
Y %T-V 


(H42) 


P 


4g 

3c 


4 ^ 

Y %T 


(H43) 


W — 

3c 


/SA A / AA. .4^ A. AA. ^ A A \ 

T TV + VT + 4y VVV) 


(H44) 


The relationships obtained for a nonisotropic photon gas are obvious general- 
izations of the isotropic case through the generalization of the scalar temperature 
concept to one of a second-order tensor. Conceptually, this generalization is 
immediately obvious . The scalar temperature is a statistical weighting parameter 
for the "conserved" quantity energy. A complex system of a large number of par- 
ticles has many more degrees of freedom than are accounted for by the single- 
particle degrees of energy , linear momentum , angular momentum , and mass . 
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TABLE H-I.- ELLIPTIC INTEGRAL $ 


a, b 




’^11^'^33 

^0 

0 ^ 0 

In $0 

1 

1 

1.00000 

1.00 

0.0000 

1 

2 

1.08904 

1.00 

.0853 

2 

2 

1.1087 

1.01 

.1032 

2 

4 

1.3139 

1.005 

.2730 

2 

8 

1.70088 

.962 

.53115 

2 

20 

2.55385 

.867 

.9376 

1 

4 

1.3440 

.978 

.2957 

1 

10 

1.9613 

.904 

.6736 

2 

10 

1.86882 

.936 

.6253 

4 

4 

1.61757 

1.020 

.4809 

4 

10 

2.12456 

.963 

.7536 

10 

10 

2.9469 

.916 

1.8675 

5 

10 

2.2728 

— 

.8210 

10 

10 

2.9469 

— 

1.8075 

10 

20 

3.9966 

— 

1:3854 

10 

40 

5.5388 

— 

1.7118 

10 

100 

8.6490 

-- 

2.157 


\)efined in equation (H16). 
= $ when V/c = 0. 




APPENDIX I 

ISOTROPIC PHOTON GAS WITH LINEAR MOMENTUM AND 
INTRINSIC ANGULAR MOMENTUM 


In this appendix , terms for an isotropic photon gas with linear momentum 
and intrinsic angular momentum are defined . 


Photon number density: 


20 0 / 4 4\ 

" = 27 K T ) 


(ID 


Total energy density: 


T c 


V ^ V ^ 

4 ^ 4 ^ 4 / 6 ''+ ^ 6 ''- 
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_ 2a_4 r. / 6 6\ 


4 4 

Y_ 


( 12 ) 


Momentum density; 


- 8 0 ( 6 " 6 ^- 
g = T ^ T + Y_ V 


3 3 
c 


(13) 


Intrinsic angular momentum density: 


0 - 20 h0 „3 

27 2, 
c k 


(v/v. 


r/v.) 


(14) 
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Total pressure: 


P - i £ 

3 c ^ 


= p5 + + g_V_ 


eV-V. 


(15) 


Local pressure: 


P = o = 


e 2aT 


3 3c 


Intrinsic angular momentum density flux: 


(y/ y v.") 


(16) 


m 


27 


(y+^ - Y_^) 6 + (2y+^ - l) 

c 


V V 


(zrj - i) 


• * - 


(17) 


Intrinsic angular momentum number density: 


m = trace (m) = |y | o T^fy/ - y_^) 


(I») 


Entropy density: 




(19) 


Entropy density flux: 


. _ 8 aT' 
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APPENDIX J 

ISOTROPIC PHOTON GAS WITH LINEAR MOMENTUM AND 
INTRINSIC ANGULAR MOMENTUM IN TERMS OF THE NET 
AND DIFFUSION VELOCITIES 


In this appendix , terms of an isotropic photon gas with linear momentum 
and intrinsic angular momentum are defined in terms of the net and diffusion 
velocities . 


Photon number density: 
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Total energy density: 
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Momentum density: 
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2 
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Intrinsic angular momentum density: 


£ = ^U + mV = |y-^g- '^ — — ■ 2 [26 + (1 + 6 u] 

^ fl- 6 ^) 


(J4) 
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Total pressure: 


P = + XT , \.Z [(l + 3J^)(w + UO) t S^(3 + S^)(VU + Cv)] 


Local pressure: 
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Intrinsic angular momentum density flux: 
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Intrinsic angular momentum number density: 


Entropy density: 
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80 IT oT^Y^a ^ 
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Entropy density flux: 
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APPENDIX K 

PHYSICAL SPACE VARIABLES FOR RANDOM POLARIZATION 


Physical space variables for random polarization are 


Photon number density: 


n 


27 kc^ 


3 4 
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Total energy density: 
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Momentum density: 
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Intrinsic angular momentum density: 


£ = 0 


Total pressure: 


P = p6 + gV 


Local pressure: 
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Intrinsic angular momentum density flux: 

m = 0 


defined in this appendix. 
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Intrinsic angular momentum number density: 


m = 0 


Entropy density: 


s 


16 a 
3 c 





Entropy density flux: 


5 = sV 


Energy-stress density four tensor: 
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APPENDIX L 

PHYSICAL SPACE VARIABLES FOR RANDOM POLARIZATION, V/c « 1 


Physical space variables for random polarization where V/c << 1 are defined 
in this appendix. 


Photon number density: 




(LI) 


Total energy density: 




(L2) 


Momentum density: 


^ _ 16 oT’ - L . 

g==-3--jV 1 + 3^ + 


(L3) 


Intrinsic angular momentum density; 


1^0 


(K4) 


Total pressure: 


P 


3 c ^ 


6(1 + 2-2 + 


V 

I — 
‘ 2 
c 


. 4 '? 1 . 2 ^ . 


~2 


(L4) 
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Local pressure: 




Intrinsic angular momentum density flux: 


m - 0 


Intrinsic angular momentum number density: 

m = 0 


(L5) 


(K7) 


(K8) 


Entropy density: 


Entropy density flux: 


s 


16 

3 


c 




(L6) 


5 = sV 


(H30) 
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